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Abstract 

If Ag{P,a,k) is the scattering amplitude, corresponding to a potential q G 
LP'{D), where D C is a bounded domain, and e^^"''^ is the incident plane wave, 
then we call the radiation pattern the function A((3) := Aq{P,a,k), where the 
unit vector a, the incident direction, is fixed, and A; > 0, the wavenumber, is 
fixed. It is shown that any function G L^(S'^), where S'^ is the unit sphere 
in M^, can be approximated with any desired accuracy by a radiation pattern: 
— A{P)\\]^2(^s2'j < e, where e > is an arbitrary small fixed number. The 
potential q, corresponding to A{P), depends on / and e, and can be calculated 
analytically. There is a one-to-one correspondence between the above potential and 
the density of the number of small acoustically soft particles Dm C D, 1 < m < M, 
distributed in an a priori given bounded domain C M^. The geometrical shape 
of a small particle Dm is arbitrary, the boundary Sm of Dm is Lipschitz uniformly 
with respect to m. The wave number k and the direction a of the incident upon 
D plane wave are fixed. 

It is shown that a suitable distribution of the above particles in D can produce 
the scattering amplitude A{a',a), a', a G 5"^, at a fixed A; > 0, arbitrarily close in 
the norm of L?'{S'^ x S"^) to an arbitrary given scattering amplitude /(a', a), corre- 
sponding to a real- valued potential q G L^{D), i.e., corresponding to an arbitrary 
refraction coefficient in D. 



1 Introduction 

Let D C be a bounded connected domain with Lipschitz boundary S. 
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The scattering of an acoustic plane wave Uq — Uo{x) — e*'^"^'^, incident upon D, is 
described by the problem: 

[V^ + A;^no(x)]M = in R^ (1) 



u = uo{x) + V, (2) 

V — A(a', a) 1- o - I , r := Ixl — > oo, — :— a'. (3) 

r \r J r 

The coefficient A{a' , a) is called the scattering amplitude, > is the wave number, 
which is assumed fixed throughout the paper, and the dependence of A on A; is not shown 

by this reason, a G iS^ is the direction of the incident plane wave, a' is the direction of the 
scattered wave, nQ{x) is the known refraction coefficient in D, nQ{x) = 1 in ]R^\D, 
and V is the scattered field. 

Let Dm-i 1 < m < M, be a small particle, i.e., 

kna <^ 1, where a—- max diamD^, kn — A; max |no(x)|. (4) 

2 l<m<M xeD 

The geometrical shape of is arbitrary. We assume that Dm is a Lipschitz domain 
uniformly with respect to m. This is a technical assumption which can be relaxed. It 
allows one to use the properties of the electrostatic potentials. Denote 

d :— mm dist{Dm, Dj). (5) 

Assume that 

a < d. (6) 

We do not assume that d ^ Aq, that is, that the distance between the particles is much 
larger than the wavelength. Under our assumptions, it is possible that there are many 
small particles on the distances of the order of the wavelength. 
The particles are assumed soft, i.e., 

u\sm = 1 < m < M. (7) 

As a result of the distribution of many small particles in D, one obtains a new material, 
which we want to be a "smart" material, that is, a material which has some desired 
properties. Specifically, we want this material to scatter the incident plane wave according 
to an a priori given desired radiation pattern. Is this possible? If yes, how does one 
distribute the small particles in order to create such a material? 

We study this problem and solve the following two problems, which can be considered 
as problems of nanotechnology. 

The first problem is: 
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Given an arbitrary function f{(3) G L^(S'^), can one distribute small particles in D 
so that the resulting medium generates the radiation pattern A{(3) := A{(3,a), at a fixed 
k > and a fixed a E S"^ , such that 



where e > is an arbitrary small fixed number? 

The answer is yes, and we give an algorithm for calculating such a distribution. This 
distribution is not uniquely defined by the function and the number 5 > 0. 

The second problem is: 
Given a scattering amplitude f{a',a), corresponding to some refraction coefficient 
n{x) in a bounded domain D, can one distribute small particles in D so that the resulting 
medium generates the scattering amplitude A{a' , a) such that 



where e > is an arbitrary small fixed number? 

The answer is yes, and we give an algorithm for calculating the density of the desired 
distribution of small particles given f{a', a), Va', a G S"^, k > being fixed. 

To our knowledge the above two problems have not been studied in the literature. 
Our solution to these problems is based on some new results concerning the properties 
of the scattering amplitudes, on our earlier results on wave scattering by small bodies of 
arbitrary shapes (see [S]), and on our solution of the 3D inverse Schrodinger scattering 
problem with fixed-energy data [7], jH], [HI- 

In Section |21 we derive some new approximation properties of the scattering am- 
plitudes. Essentially, we prove the existence of a potential q G L'^{D) such that the 
corresponding to this q scattering amplitude v4(/3), (3 = a', at an arbitrary fixed a G 5^ 
and an arbitrary fixed k > 0, approximates with any desired accuracy any given function 
/(/?) G L?{S'^) with a small norm. Moreover, we give formulas for calculating this g, 
and these formulas often work numerically for / which are not small. The potential q is 
related explicitly to a certain distribution of small particles in D. Consequently, we give 
formulas for calculating this distribution. 

In Section El we drive an equation describing the self-consistent field in the medium 
consisting of the small particles distributed in D. This equation is equivalent to a 
Schrodinger equation with a potential q{x) supported in the bounded domain D and 
related in a simple way to the density of the distribution of the small particles. 

The author has solved the 3D inverse scattering problem of finding a compactly 
supported potential q from the knowledge of noisy fixed-energy scattering amplitude [H]. 
This algorithm allows one to calculate qsi^x) from the knowledge of noisy data /^(q;',^), 
suPa',ae52 1/5 - /I < such that 



||/(/?)-A(/?)lU.(5.)<£, 
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||/(a',a) - A(a',a)||L2(52x52) < e 



(9) 



sup \qs{x) - q{x)\ < ri{6) 








(10) 



3 



where q{x) is the exact potential, generating the exact scattering amplitude f{a', a) at 
a fixed A; > 0. 

Applying this algorithm to the exact data /(a', a) or to the noisy data fs{a', a), one 
obtains a stable (in the sense (fTUj) ) approximation of q, and, consequently, of the density 
of the distribution of small particles, which generates the scattering amplitude arbitrarily 
close to the a priori given scattering amplitude. 

The author's solution of the 3D inverse scattering problem with the error estimates 
is presented in Section HJ 



2 Approximation properties of the scattering ampli- 
tudes 

If A; > is fixed, then the scattering problem (HJ-© is equivalent to the Schrodinger 
scattering problem on the potential qo{x): 

[V^ + k'^ ~ qoix)]u = inM^ (11) 



in D', D' ■.= R^\D 
k'^[no{x) — 1] in D. 



'12) 



The scattering solution u = u^^ solves (uniquely) the equation 

r ^ik\x-y\ 

Uqo = uo- g{x,y)qo{y)uq,Xy)dy, 9{x,y) := — r 

Jd 4:7^\x-y\ 



The corresponding scattering amplitude is: 



Ao{a',a) = -— / e 



-ika'-x 



qo{x)uqo{x,a)dx, 



where the dependence on k is dropped since > is fixed. 

If go is known, then Aq := Ag^ is known. Let q G L'^{D) be a potential and Aq[ 
be the corresponding scattering amplitude. Fix a E S"^ and denote 



A{j3) := Aq{a',a), a' = p. 



Then 



A{l3) = --^ [ e ''''^'''h{x)dx, h{x) := q{x)uq{x,a). 
Theorem 1. Let f{(3) G L'^{S'^) he arbitrary. Then 



inf 

h&L'^{D) 



An 



''''^■''h{x)dx 



D 



0. 



(13) 

(14) 

a', a) 
(15) 

(16) 
(17) 
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Proof of Theorem{J\ If fITTj) fails, then there is a function G L'^{S'^), / 7^ 0, such 
that 

/" /" e-^'=^-^/i(a:)dx = WheL^D). (18) 

This imphes 

ip{x) := [ dl3f{l3)e-'^^-'' = Vx e D. (19) 
The function v?(x) is an entire function of x. Therefore (|19|) implies 

ip{x) = Vx G Ml (20) 

This and the injectivity of the Fourier transform imply f{j3) = 0. Note that f{x) is the 
Fourier transform of the distribution f{f])5{k — X)X~'^, where S{k — X) is the delta-function 
and Xf5 is the Fourier transform variable. The injectivity of the Fourier transform implies 
f{P)5{k - A) = 0, so fiP) = 0. 

Theorem 1 is proved. □ 

Let us give an algorithm for calculating h{x) in (jl7|) such that the left-hand side of 
(fTTj) does not exceed e, where e > is an arbitrary small given number. 

Let {y^(/5)}^0' ~ ^,m, —i<m<i,he the orthonormal in L'^{S'^) spherical 
harmonics, 

Ye,U-P) = = {-lY^"'Ye,Ul3), (21) 

where Jf are the Bessel functions. It is known that 



e=o-e<m<e ' ' 

Let us expand / into the Fourier series with respect to spherical harmonics: 

f{P)= E kmYUP). (24) 

e=o,-e<m<E 

Choose L such that 

J2\fi,n\'<e'. (25) 

e>L 

With so fixed L, take hi^mij), < £ < L, —i < m < i, such that 

7r\i/2 



= -(-^)' (^) J r'/^Je^{kr)he,m{r)dr, (26) 
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where 6 > 0, the origin O is inside D, the ball centered at the origin and of radius h 
belongs to D, and hs^rair) — r > b. There are many choices of h^rn{f) which satisfy 
ira . If (123) and (I2S1) hold, then the norm on the left-hand side of ^ is < e. 
A possible explicit choice of he^rn{r) is 

e<L, 



ti,m=s (-^)Vife^M+^W^ (27) 
0, i>L, 

where we have assumed that 6 = 1 in ()26|). and used the following formula (see formula 
8.5.8]): 

1~ 



f28) 



9^^Ak), 



where S^^k) are Lommel's functions, r(a;) is the Gamma-function, hi^rn{r) in ()27|) do 
not depend on r, and we assume that h{x) = for r := |x| > 1. 

Let us prove that for any q & LF' there exists a g G L^{D) such that q{x)uq approxi- 
mates h{x) in L^(D)-norm with arbitrary accuracy. 

Theorem 2. Let h G L^{D) be arbitrary. Then 

inf \\h-qu,{x,a)\\=0. (29) 

Here a E and k > are arbitrary, fixed. There exists a potential q such that h = qu, 
^/ 1 1^1 is sufficiently small. 

Proof of Theorem{^ In this proof we first assume that the norm of / is small, and then 
we drop the "smallness" assumption. If the norm of / is sufficiently small, the the norm 
of h is small, so that the condition 

inf |Mo(a;)- / g{x,y)h{y)dy\> (30) 

x€D Jjj 

is satisfied. Here g is defined in formula (jl3|) . If this condition is satisfied, then the 
formula 

q{x) = h{x)[uo{x) - g{x,y)h{y)dy]-^ (31) 
Jd 

yields the desired potential q. The function h generates the function u := ^, where u 
is the scattering solution, corresponding to the potential q, constructed by formula 
Therefore, the infimum in ()29p is attained if condition (jHUj) is satisfied by the given h. 

If / is arbitrary, not necessarily samll, then h is not necessarily small. If, nevertheless, 
condition (IHUI) holds for this h, then the potential q, given by formula ((SH), belongs to 
L'^{D) and yields the scattering amplitude Ag{P) which satisfies ((HI). 
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On the other hand, if condition (j3(J|) does not hold, then formula (|3ip may yield a 
potential which is not locally integrable. In this case, as was proved in ^U], one can 
perturb h slightly, so that the perturbed h, denoted by h^, \ \h — he\\L2{D) < would yield 
by formula ()31|) . with he in place of h, a potential Qs G L'^{D). This potential generates 
the scattering amplitude Aq^{P) which satisfies estimate (jH)), possibly with ce in place of 
e, where the positive constant c does not depend on e. Theorem 2 is proved. □ 

Let us give a different point of view on the role of the "smallness" assumption. If 
||'Z||l2{d) 0, then the set of functions qu becomes a linear set. Thus, if fails, then 
there exists an /i 7^ 0, /i G L'^{D) such that 

h{x)q{x)ug{x,a)dx = \/qeL^{D) (32) 

JD 

Condition (j32j) holds in the limit H^H — > because in this limit the set of functions qu 
becomes linear. 

Let c = const > be small. We will take c — > eventually. Choose 

For sufficiently small c > the equation 

Jd 

is uniquely solvable for Ug in C{D) because ||T|| < 1 if c > is sufficiently small. We 
have 

quq = qe''"'-^ - qTuq = ch + O(c^), c ^ 0. (33) 
Substitute into divide by c, and take c — * 0. The result is: 

\h\'^dx = 0. (34) 



D 



This implies h = 0. 

We describe the relation between q{x) and the density distribution of small particles 
in Sectional This relation makes it clear that a suitable distribution of small particles will 
produce any desirable potential q G L'^{D), and, consequently, any desirable scattering 
amplitude (radiation pattern) at an arbitrary fixed a G S*^ and A; > 0. 

We describe the algorithm for calculating the above distribution of small particles, 
given /(/3) G L\S'^), in Section El 

3 Scattering by many small particles. 

The exact statement of the problem is: 

M 

[V' + e- qo{x)]u = in R3 \ y Dm, (35) 



m—l 
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on IJ Sm, Sm ■■= dDm- (36) 



M 

u 

m=l 

u = e*'^"-^ + v:=Uo + v, (37) 



e*'"' / i \ , , , X 







ay 


r := \a\ 







r 



(3J 



We look for the solution of the form 



M 

u{x) = uq{x) + / G{x, s)am{s)ds, (39) 

m=l '^•5''" 

where G{x, s) is the Green function which solves the scattering problem in the absence 
of small particles, i.e.: 

[V^ + A;^ - qo{x)]G = -6{x - y) in M^ (40) 

lim \x\(^- tkG] = 0, (41) 
I^Hoo \d\x\ J 

and Uq is the corresponding scattering solution. It was proved in [3], [7], that 

^ik\x\ 



'^(^^y) = :r~n'^o{y, a) + { —r ] ,\x\ ^ oo, a = -—, (42) 



where Uq is the scattering solution corresponding to go- 

The function solves equation (jHHjl . satisfies the radiation condition (jHH|) . because 

^^0 = ^^0 + ^^0, (43) 
where Vq satisfies the radiation condition (PT|) . Therefore (jH^ solves the problem ()35|)- 



if dm are such that the boundary condition ()36|) is satisfied. All the above did not 
use the smallness of the particles. 

Let us now use assumptions (@)) and (jH)). Let Xj G Dj be an arbitrary point inside 



Dj. Then 



sup \G{x, s) — G{x, Xj) \= O (ka + -j , \x — Xj\ > d. (44) 



This follows from the integral equation, relating G and g: 

G{x,y) = gix,y) - / g{x,z)qo{z)G{z,y)dz, 



8 



and from the estimates: 



gifc|x— s| 




1 3^ *^ 1 




1 


^ik(\x—s\~ 


-\X-Xj\) 


1 J 1 


An 


X — s 


An 






An 








X — s\ 





k 



k\x-xj\ (l + O +0(A;a; 



X — s\ 



\x — Xj — {s — Xj)\ = \x — Xj\ ( 1 + O ( - 



From the integral equation for G it follows that 

Gix,y) ^ g{x,y) as x ^ y. 
Therefore one may approximate (|39p as 



M 



u(x) = Uo{x) + G{x, Xj)Qr, 



m=l 

where \x — Xj^ \ > d for all m, 1 < m < M, and 



l + 0\ ka+ - 
a 



Qr 



arn{s)ds. 



(45) 



(46) 



Therefore, if one knows the numbers Qm-, 1 < m < M, then one knows the scattering 
solution u{x) at any point which is at a distance > d from the nearest to x small body. 

Let us derive a linear algebraic system for calculating Qm- To do this, let us use the 
boundary condition (jHBjl . We have: 



/ G{s,t)amdt = - 



Since k\s — t| < 2a ^ 1, one has 



G(s,t) 



^ik\s—t\ 



An\s — t\ An\s — t 
Consequently, equation (jlTjl can be written as 

o"m(i) 



Uo{Xm) + ^ G{Xm, Xj)Qj 



'l + 0{ka)). 



(47) 



47r|s - t\ 



-dt 



^o('^m) ~l" ^ ^ C('^m) •^j)Qj 



(4J 



(49) 



This is an equation for the electrostatic charge distribution am on the surface Sm of the 
perfect conductor Dm, charged to the potential which is given by the right-hand side of 
(^nj- The total charge on the surface of the conductor is given by the formula: 



)Q, 



1 < m < M, 



(50) 
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where Cm is the electrical capacitance of the conductor D^- Equation (|5Up is a linear 
algebraic system for the unknown Qj. 

Assume that the distribution of small bodies Dm in D is such that 



lim > Cm = / C(x)dx, 



M 

D„,CD 



(51) 



where D is an arbitrary subdomain of D. This means that C{x) is the limiting density 
of the capacitance per unit volume around an arbitrary point x G -D. Then the relation 
fjIKjl in the limit 

a 

takes the form 

u{x) = uo{x)- G{x,y)C{y)u{y)dy, (52) 
Jd 

where C{x) is defined in (jS^- An equation which is similar to ()52p . with g{xj, y) in place 
of G{xj,y), has been derived in 2J by a different argument. 
Equation is equivalent to the Schrodinger equation 

[V^ + e- qo{x) - C{x)] u = 0, (53) 

and u{x) is the scattering solution corresponding to the potential 

q{x) = qo{x) + C{x). (54) 

If ^'0(3;) is known (which we assume), then q{x) and C{x) are in one-to-one correspon- 
dence. 

If the small particles Dm are identical, and Co is the electrical capacitance of a single 
particle, then 

C{x) = N{x)Co, (55) 

where N{x) is the density of the number of particles in a neighborhood of the point x, 
that is, the number of particles per unit volume around point x. 

Therefore, given /(/3) G L'^{S), one finds q{x), such that ||Aq(/5) — f{/3)\\ < e, where 
Ag{(3) is the scattering amplitude, corresponding to the potential q, the energy k"^ > 
and the incident direction a being fixed, and /3 = a' is the direction of the scattered 
wave. 

Let us describe the steps of our algorithm. 
Step 1. Given /(/5), find h G L^iD). 

Step 2. Given h G L'^{D), find q such that \\h — q{x)uq{x)\\L2(^£)^ < e. 

Let us elaborate on Step 2. First, assume the existence of a potential q, such that 
h = qu. Consider the equation 

u = uo- j gquqdy = ^0 - j ghdy. (56) 
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We have 

qUq := h. 

Thus, 



D 



Ait 

Multiply by q. Then 

h = Uoq-q I ghdy 

Therefore, if 



Ag{P) = I e-'^^-''h{x)dx. (57) 



D 



inf \uo{x) - I g{x, s)h{y)dy\ > 0, 

x&D 

then the solution of the equation qug = h is unique and is given by the formula: 

/ N h(x) , , 

"0(2;) - J D ghdy 

Formula (|58|) yields a potential for which Aq{P) is given by formula (j57jl . and the cor- 
responding scattering solution is given by formula ()5fij) . All this is true provided, for 
example, that 



sup 



g{x,y)h{y)dy 



D 



< 1. (59) 



Inequality ()59|) holds if h is fixed and diam D is sufficiently small, because of the following 
estimate: 



sup| / g{x,y)h{y)dy\ < (Ait) 2\\h\\L2(^r)){diam D)^ . 



D 



Inequality (j59|) also holds if ||/i||L2(£)) is sufficiently small and D is fixed. The norm 
||^||l2(d) is small if ||/||i2(52) is sufficiently small. For the formula to yield the 
desired potential, the inequality is not necessary. If one can find a potential q{x) 
from the given h by formula ()58|) . then this q generates the scattering solution by the 
formula 

Uq = uo- / ghdy (60) 



D 



and 

h = q{x)uq{x). (61) 

The potential q can be found by formula (j3Hj) . provided that /(/?) is sufficiently small, 
because then h will be sufficiently small as follows, e.g., from ()27|). 
If q is found, then 

Nix) = ^MZ^. (62) 

Thus, the corresponding distribution density of small particles is given analytically. 

Analytical formulas, which allow one to calculate Cq with any desired accuracy, are 
derived in see also formula (91) below. 
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Remark. If f{j3) corresponds to a real-valued q{x) , then formula (|58|) yields a real-valued 
potential. In general, formula (j3S|) yields a complex-valued potential. To get a complex- 
valued potential by a formula, similar to ()55|) . one has to replace the boundary condition 
(jHUI) by the condition 

un = (u on Sm, (63) 

where N is the exterior unit normal to the boundary S, and C, is a complex constant, the 
impedance. Then Cq in (f^ should be replaced by the quantity: 



c 



= (64) 



(see fS^), and, therefore, formula (jKHjl yields a complex-valued C(^{x) if ( is a complex 
number. 

Suppose that a given h corresponds to a potential q{x) G L'^{D) in the sense that 
h = q{x)u{x), where u{x) is the scattering solution corresponding to this q{x) at the 
wavenumber k > and with the incident direction a. Then formula (58) defines q{x), 
and the corresponding scattering solution is -u = 

If formula (58) does not produce a p G L'^{D), then one can replace h in (58) by an 
hir, \\h — /i£||L2(-£)) < e, and get a square-integrable potential q^ by formula (58) with h 
replaced by h^. If e is sufficiently small, this potential q^ generates the radiation pattern, 
which differs by 0{e) from the desired /. 

4 Ramm's solution of the 3D inverse scattering prob- 
lem with fixed-energy data 

We follow [0], [ZI. Consider first the inversion of the exact data Aq{a',a). Let 

oo 

Ag{a',a)= Ai^rn,e',m'Yi'^rn'ia')Yi^mia)- (65) 

£,£'=0 



1 „r\Im8\ 

\Y,{0)\<^-^, Vr>0, 9eM, (66) 



It is proved in [7j that 



where 

3 

M = {z : z e C^,z ■ z = k'^}, z ■ ( ■.= ^ ZjQ. 

i=i 

Estimate (66) allows one to prove ( P) that the series (65) converges absolutely for 
a' = 9' E M., so that the exact data Ag{a', a) allow one to calculate the values Ag{6', a), 
9' E M.. These values are used below in the inversion formula (68). 
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One can prove [7j, that any ^ G can be written (nonuniquely) as 

^ = e'-e, e',9eM, i^i^oo. (67) 

In jZI explicit analytical formulas are given for 6' and 9 satisfying ()67|1 . 

The exact data A{a',a) admit an analytic continuation from 5*^ x 5*^ onto Ai x S'^. 

Let 

Jd 

The inversion formula, proved in [7j, is 

^ ' ^ (68) 



g(^) = lim 



-An / A{e' ,a)i'{a,e)da 



where ffHTj) holds and z^(q;, ^) is an arbitrary approximate solution to the problem 

J^(z/) := / \p{x)\'^dx = inf ■.= d{e). (69) 

Jai<\x\<b 



' ai<|a;|f 

Here 

p[x) ■= e~'^ '' I u{x,a)u{a,6)da - 1, (70) 

Oi > is a radius of a ball which contains D as a strictly inner-subdomain, and b > ai 
is an arbitrary fixed number. The approximate solution ly to (|69|) is understood in the 
following sense: 

< 2d{e). (71) 

This means that it is not necessary to find a very accurate approximation of the infimum 
in problem (69). It is sufficient, for example, to find any function i/(a,^) for which the 
functional (69) takes the value not more than 2d{6), and the inversion formula (68) holds 
with such i>{a,9). Also, formula (73) below, with the error term, holds as well. 
It is proved in jTj that 

d{9)<^^, 9eM, (72) 

where c = c(||g||) > is a constant depending on an L°°{D) norm of q. Therefore, given 
the exact data Ag{a',a), one recovers the potential q{x) by formula (jUHj) . 
The error estimate of formula ()68p is given by the formula: 

(73) 



q{0 = -47r Ai9', a)z/(«, 9)da + O (^-^^ , \9\ ^ oo, 



where (jU7j) holds. 

If q{x) is found, then 

N{x)Co = qix)-qo{x), (74) 
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so that the density of distributions of small particles is found explicitly. 
Consider now inversion of noisy data As{a', a), 

sup \As{a' , a) — A{a' ,a)\ < 6. 

a' ,a(^S'^ 



(75) 



Here A{a\ a) corresponds to an exact potential, and A{a\ a) is not known. Instead, its 
noisy measurements As{a',a) are known. 
Define 

|ln5| 



N{5)- 

where [x\ is the integer nearest to x > 0, 



_ln|ln(5| 



N(S) 



l,(^',a) = ^A,,(a)F,(^'), E = E E ' 



£=0 

N(5) 



£ -l<m<l 



usix, a) = e^'^"-^ + Asiia)Yi{a')heikr), a' := -, 



r = \x\ 



£=0 



ps[x; V) = e 



us{x,a)u{a)da — \, 9 E A4, 



52 

7Af(<5) 



/i(<5) = e~ 
bo '■= 2 d'i-o,'>TT- D, 



7 = ln-i >0, 
bo 



K = \Im9\ 



Let 



bo < ai < b, 

where ai and b are arbitrary positive fixed numbers. Consider the problem: 

1^1 = sup := i}{6) 

under the constraints 

1^1 [\\ps{^)\\Lmx:a^<\x\<b}) + \W\\ L^S^)^''^ M] < C, 



(76) 



(77) 

(78) 

(79) 
(80) 



f82) 



(83) 



OeM, 9' -6 = ^, 9',9eM, (85) 



where c > is a sufficiently large constant, and 60 < < b. 
It is proved in [7] that 
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Let 6{5) and ^'^(a) be any approximate solution to (jHHj) -(|S ^ in the sense that 

m\>^m- (87) 

Define 

qg-=-4:n As{e',a)us{a)da. (88) 

The following result is proved in 
Theorem (Ramm). One has 

sup|*-fK)l = o(^i!^). (89) 

This result gives an inversion formula for finding the potential from noisy fixed-energy 
scattering data. 

Thus, the algorithm for finding the density of the distribution of small particles from 
the fixed-energy scattering data A{a', a) can be formulated as follows: 

Step 1 . Given A{a',a), find q{x) using the inversion formulas in the case of the 
exact data or (jHHj) in the case of noisy data. 

Step 2 . Find the density of the distribution of the small particles by formula (jHS)), where 
formulas for Co are given in [7]. 

|Co-C(")| = 0(2"), 0<Q<1, (90) 

where Q depends only on the geometry of the surface, 

C^"^=4^l^l'(feT^ / / ^ / / ^P{t,tl)...Htn.,Qdt^■dt^ (91) 

{ [Zir)"- Js Js Tst J g n integrals J g J 

d 1 

ilj{t,s) = —— — , rst = \s-t\, \S\ = measS, (92) 

S is the surface of the conductor, Cq is the electrical capacitance of this conductor, and 
Nt is the exterior normal to S at the point t. 
In particular, for n = one gets 

cm = 1^, J..^n^J±. (93) 

J J S Js ^st 

It is proved in [S] that 

C(°) < Co. (94) 

Formula (PT|) given an approximate value C*^"-* of the electrical capacitance of a perfect 
conductor placed in the space with dielectric permittivity ^o = 1- If 7^ 1; then one has 
to multiply the right-hand side of (PT|) by Sq. 
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